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Abstract 

A polynomial complexity algorithm is designed which tests whether 
a point belongs to a given tropical linear variety. 


Introduction 

Consider a linear system 

A-X = b (1) 

with the m x n matrix A = (a^j) and the vector b = (5*) defined over the 
field K = C((f 1//o °)) = {c = c 0 t l °^ q + c 1 t < ' l0+1 ^ q + • • •} of Puiseux series where 
i 0 G Z, 1 < q G Z. Consider the map Trop(c ) = i 0 /q G Q and Trop( 0) = oo. 
Denote by P C K n the linear plane determined by the system (JT]) . The closure 
in the euclidean topology Trop(P) C M” is called a tropical linear variety [13] 
(for the basic concepts of the tropical geometry see [12], jl~3]). 

More generally, the tropical variety attached to an ideal / C K[X U ..., X n \ 
is defined as Trop{U) C W 1 where U C K n is the variety determined by 
I. A tropical basis of / is a finite set fi, ■ ■ ■, fk £ I such that Trop{U ) = 
Trop(V(fi )) fl • • • fl Trop{y{fk)) where V(f \) C K n denotes the variety of all 
the zeroes of f±. In [3] an algorithm is devised which produces a tropical basis 
of an ideal. Having a tropical basis, one can easily test, whether a point v 6 M" 
belongs to the tropical variety Trop{U ) since Trop{V{f )) = V(Trop(f )) (due 
to ra) where the tropicalization Trop(f) of a polynomial / G K[X- l, ..., A n ] 
is defined coefficientwise, and V(Trop(f)) is the tropical hypersurface of all 
the zeroes of the tropical polynomial Trop(f). But on the other hand, in 
[3] an example is exhibited of a tropical linear variety with any its tropical 
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basis having at least an exponential number of elements, while the algorithm 
recognizing a tropical linear variety designed in the present paper has the 
polynomial complexity. 

We study the problem of the complexity of recognizing Trop(P). In other 
words, we design a polynomial complexity algorithm which for a given vector 
v = (rq,... ,v n ) 6 (If! Q) n with real algebraic coordinates tests, whether a 
system ([I]) has a solution x — (x±,..., x n ) G K n with Trop(x ) = v. Obviously, 
this captures also the case v G ((R flQ) U {oo}) n , so we can w.l.o.g. assume 
below that v G (R D Q) n . 

Observe that the problem of recognizing, whether just the zero vector be¬ 
longs to a given tropical (non-linear) variety, is ATP-hard (see |SJ) since the 
solvability of a system of polynomial equations from Q[Xi,... ,X n ] is equiva¬ 
lent to that the zero vector belongs to the tropical variety determined by this 
system. 

We mention also that testing emptiness of a tropical non-linear prevariety 
(i. e. an intersection of a few tropical hypersurfaces) is NP- complete [TS], while 
testing the emptiness of a tropical linear prevariety belongs to NP D coNP 

(H. Bl. HU)- 

1 Algorithm lifting a point to a tropical linear 
variety 

We assume that the entries a t j, b t G K, 1 < i < m, 1 < j < n are provided 
in the following way (cf. [5], [9], [2j). A primitive element z G K is given as 
a root of a polynomial equation h(t,z) = 0 where h G Z[£, Z\, and by means 
of a further specifying a beginning of the expansion of 0 as a Puiseux series 
over the held Q of algebraic numbers (to make a root of h to be unique with 
this beginning of the expansion). Also we are supplied with rational functions 
hi : j,hi G Q(t)[Z] such that cqj = \j(z), bi = hj(z), 1 < i < m, 1 < j < n. 
We suppose that deg(/r), deg(hjj), deg(hj) < d. In addition, we assume that 
each rational coefficient of the polynomials h, hij, hi is given as a quotient of a 
pair of integers with absolute values less than 2 A/ . The latter means that the 
bit-size of this rational number is bounded by 2 M. 

First, the algorithm cleans the denominator in the exponents of the Puiseux 
series of z replacing t x ' q by t for a suitable q < d to make z to be a Laurent 
series with integer exponents (and keeping the same notation for z, h, h^j , hi). 
Here and below to develop a Puiseux series with coefficients from Q within 
the polynomial (in the bit-size of the input and in the number of terms of 
the expansion) complexity, we exploit the algorithm from [6]. The algorithm 
makes use of a presentation of the held of the coefficients as a finite extension 
of Q via its primitive element (see e. g. 0; 05 0) similar to the presentation 
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of o.j j , b, above (we don’t dwell here on the details since it does not influence 
the main body of the algorithm, for the sake of simplifying the exposition a 
reader can suppose that the coefficients of the Puiseux series of z are rational). 

The coordinates of the vector v we also multiply by s and keep the same 
notation for v. We say that two coordinates v n , v 32 of v are congruent if 
v ji ~ v n £ Every class of congruence provides a subsystem of (JT]) (note that 
just one of such classes can contain the column b, namely the class of Vj being 
integers). Observe that these subsystems are pairwise disjoint (supported in 
disjoint columns), and their conjunction is equivalent to ([T]). 

We assume that the vector v is provided in the following way (cf. |10j . ||2] 
and also above). A primitive real algebraic element u G Q D M is given as 
a root of a polynomial g G Z[Y] together with specifying a rational interval 
[ei,e 2 ] which contains the unique root u of g. In addition, certain polyno¬ 
mials gj G Q[Y], 1 < j < n are given such that Vj = gj(u). We suppose 
that deg(g), deg(c/j) < d and that the absolute values of the numerators and 
denominators of the (rational) coefficients of g , g 3 and of ei,e 2 do not exceed 
2 m . 

To detect whether for a pair of the coordinates the congruence Vj l — v 32 G Z 
holds, the algorithm computes an integer approximation e G Z of \vj 1 —Vj 2 —e\ < 
1/2 (provided that it does exist) with the help of e. g. the algorithm from [2] 
and then verifies whether v n —Vj 2 = e exploiting 0.0 or pj. This supplies us 
with the partition of the coordinates v±,, v n into the classes of congruence. 
Thus, fixing for the time being the subsystem of flU) corresponding to a class 
of congruence, one can assume w.l.o.g. that the coordinates of v are integers. 
Moreover, one can assume w.l.o.g. that v = 0 replacing each a h j by t~ Vj ■ a i)3 
(and keeping for them the same notations). 

Then by elementary transformations with the rows of matrix A over the 
quotient-ring Q (t)[Z]/h (making use of the basic algorithms e. g. from [2j) and 
an appropriate permutation of columns, the algorithm brings A to the form 
a,i t i = 1, a^j = 0, 1 < i j < m (one can assume w.l.o.g. that rk(A) = m). 

For m < j < n denote r 3 := — min i<i< m {Trop(a it j)}. If r 3 < 0 we put the 
coordinate x 3 — 1. Else if r 3 > 0 we put x 3 = y 3 y + y 3} \ - t- 1— • + y 3 ^ :i ■ t rj with 
the indeterminates y 3 y ,..., y^ rj over Q. 

For the time being fix 1 < i < m. Denote s* = 

min m< j< n {Trop(aij), Trop{bi)}. The i-th equation of (JTJ) one can rewrite as 
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For every Sj < k < 0 one can express the coefficient of J2 m<3<n cii,j ■ x 3 — bi 
at the power t k as a linear function L^ over Q in the indeterminates Y := 
{Vj,h m < j < n, 0 < l < r 3 }. 
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Consider the linear system 

L i} k = 0, 1 < i < m, Si < k < 0 (3) 

in the indeterminates Y. The algorithm solves (J3]) and tests whether each of 
n linear functions from the family 

L \= {L i0 , 1 <i <m; y j>0 , m< j < n} 

does not vanish identically on the space of solutions of ([3]). If all of them do 
not vanish identically then take any values of Y which fulfil (J3J) with non-zero 
values of all the linear functions from the family L. Then the equation f[2l) 
determines uniquely Xi with Trop(xi ) = 0(= vf) and provides a solution x 
of the system (JTJ) satisfying Trop(x ) = v. Otherwise, if some of the linear 
functions from the family L vanishes identically on the space of solutions of 
(J3J) then the system (JT() has no solutions with Trop(x) = v. 

To test the above condition of identically non-vanishing of the linear func¬ 
tions from the family L the algorithm finds a basis W\, ... ,w r G (Q) N and a 
vector w G (Q) N where N — |Y| such that the r-dimensional space of solutions 
of 03]) is the linear hull of the vectors W\, ... ,w r shifted by the vector w. If 
each linear function from the family L does not vanish identically on this space 
then all of them do not vanish on at least one of the vectors from the family 

F \= {w + p l ■ wi, 1 < p < nr + 1} 

1<1<R 

because any linear function can vanish on at most of r vectors from F due to 
the non-singularity of the Vandermond matrices. So, the algorithm substitutes 
each of the vectors of F into the linear functions from L and either finds a 
required one Y or discovers that ([1]) has no solution with Trop(x) = v. 

To estimate the complexity of the designed algorithm observe that it solves 
the linear system ([3]) of the size bounded by a polynomial in n, d with the 
coefficients from a finite extension of Q having the bit-size less than linear in M 
and polynomial in n, d (again for the sake of simplifying the exposition a reader 
can think just of the rational coefficients, cf. the remark above at the beginning 
of the present Section). Thus, the algorithm solves this system within the 
complexity polynomial in M, n, d (see e. g. 0). by a similar magnitude one 
can bound the complexity of the executed substitutions, and finally we can 
summarize the obtained results in the following theorem. 

Theorem 1.1 There is an algorithm which for a tropical linear variety 
V := Trop(P) defined by a linear system (QJj over the field K of Puiseux 
series, recognizes whether a given real algebraic vector uG ((KflQ) U (oo}) n 
belongs to V. If yes then the algorithm yields a solution x G K n of (Q]) with 
Trop(x ) = v. The complexity of the algorithm is polynomial in the bit-sizes of 
the system (QJ] and of the vector v. 
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Further research 


Let a point v G M n don’t lie in a tropical linear variety V = Trop(P) (cf. the 
Introduction). Theorem 11.11 implies that one can verify the latter within the 
polynomial complexity. Since hi, 0 entail that there exists a finite tropical 
basis /i, ..., fk G / of the ideal / C K[X i,..., X n \ of P. This allows one to 
find / G {/i,..., fk} C / such that v is not a tropical zero of the tropical 
polynomial Trop(f). But the number k can be exponential ([3j). Is it possible 
to construct / 0 G / for which v is not a tropical zero of Trop(fo), within the 
polynomial complexity? 

What is the complexity to detect whether a given tropical linear prevariety 
V {R \,..., R p ), with Ri,... ,Rp being tropical linear polynomials, is a tropical 
variety? 
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